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Abstract 

It is well known that any multimode positive definite quadratic Hamil- 
tonian can be transformed into a Hamiltonian of uncoupled harmonic os- 
cillators. Based on this theorem, the multimode thermal squeezed coher- 
ent states axe constructed in terms of density operators. Decoherence of 
multimode thermal squeezed coherent states is investigated via the char- 
acteristic function and it is shown that the decohered (reduced) states are 
still thermal squeezed coherent states in general. 
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1 Introduction 


In the past three decades, decoherence (reduction) has become more and more 
widely recognized, being studied in many fields of physics from non-equilibrium 
statistical mechanics [l] to quantum measurement and quantum cosmology [2, 3]. 
The purpose of this paper is to investigate decoherence of multimode thermal 
Squeezed Coherent States (SqCS’s). 

In the literature there are many equivalent definitions for one-mode, two- 
mode [4, 5, 6, 7] and multimode SqCS’s [8]. However, thermal SqCS’s are usually 
defined for one-mode [9]. (Two-mode thermal SqCS in thermo-field formalism 
is effectively one-mode.) Therefore, in this paper we first introduce a general 
definition of the multimode thermal SqCS’s in terms of density operators. Then 
we will discuss two related representations — the Wigner function and the char- 
acteristic function — and will show that the latter is the better representation for 
decoherence problems. Finally we will use the characteristic function to study 
the decoherence of multimode thermal SqCS’s. 

This paper is organized as follows: In Sec. 2 notations, conventions and a 
lemma on matrix are introduced for the mathematics used in this paper. In 
Sec. 3 a unified definition of multimode SqCS’s with the aid of a special kind of 
Hamiltonian is presented. In Sec. 4 the multimode thermal SqCS is constructed 
by thermalizing the multimode SqCS defined in Sec. 3. In Sec. 5 we calculate 
the Wigner functions and the characteristic functions of some multimode thermal 
SqCS’s. In Sec. 6 the decoherence of multimode thermal SqCS is effectuated and 
it is shown that the decohered state is still a thermal SqCS. 


2 Mathematical Preliminaries 

Throughout this paper, fi is set equal to 1; “f” denotes hermitian conjugate and 
“t” denotes the transpose of a matrix. The physical system under consideration 
is of n degrees of freedom, hence the dummy indices run from 1 to n unless 
otherwise specified. 

We use x =< x \ , X 2 , . . • , x n > and k =< Arj, fc 2 , . . . , k n > for the n-dimensional 
canonical coordinate and momentum respectively. Thus < x; k > is a vector in 
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2n~dimensional phase space, q and p denote the n-dimensional position and 
momentum operators corresponding to the canonical variables x and k. The 


Canonical Commutation Relations (CCR’s) are: 

[ft , Qj] = \Pi,Pj] = 0, [ft, Pj] = iSi 3 • ( 1 ) 

|0) denotes the n-mode Fock vacuum state, i.e., the ground state of an n- 
dimensional harmonic oscillator with unit mass and frequency: 

<x|0> = tt - * exp[-t(f) 2 ]. (2) 

The number operators axe defined in the ordinary way: 

Ni = \ tf? + «?-!)• (3) 

The (phase space) displacement operator D{< x; k >) is defined as: 

D(< x; k >) = exp[i(£- q — x • p)]. (4) 

D(< x; k >) is unitary and has the following properties: 

D\< x;k>) = b~\< x-k>) = D(- < x; £ >), (5) 

D(< x\k >) <q;p> Z) -1 (< x; k >) =< (q — x); (p — k) > . (6) 

The coherent state is defined as [4]: 

|x,£) = Z)(<x,£>)|0). (7) 


Another kind of unitary operator we will use in this paper are the elements 
of metaplectic group Mp(2n,R) — the quantum analogue of symplectic group 
Sp(2n,R). Mp(2n,R) is an n(2n + l)-dimensional Lie group with its algebra 
spanned by {qiqj,PiP 3 , qiPj +PjQi}- The elements of the Lie algebra of Mp(2n,R) 
can be organized as anti-hermitian operators in the following form: 


^ 2 

$( m ) = nE + frjPtPj + 7u(<7 iPi + p&)\ 

1 1 

i ? ? ( <* l\ ? ? t 

= 5 <9;p> U‘ p) <v ' r> 

= i < Tm < q-p >*, 


(8) 
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where = otji , (3 y — flji and 


m — 




G sp(2n,r) 


is a 2n x 2n real matrix [10], while 


( 9 ) 


J = ] , 1 = n x n unit matrix. 

V-1 OJ 

From CCR’s, we have: 

[<l(m),< q]p >*] = 


< <?;p >*= < <7;p > f i 


( 10 ) 


(ii) 


and 

[l(m 1 ),l(m 2 )] = $([m 1) m 2 ]). (12) 

Therefore the Lie algebra of Mp(2n,R) is isomorphic to sp(2n,r) — the Lie al- 
gebra of Sp(2n,R) 

The action of exp[4>(m)] E Mp(2n,R) on < q;p > can be defined and calcu- 
lated from (11): 


exp[<l(m)] < q;p>* exp[— $(m)] = exp(— m) < q\p >*, (13) 


where exp(— m) E Sp(2n,R). 

Now we replace exp(— m) in (13) by a general element S E Sp(2n,R) and try 
to find a unitary operator U(S) E Mp(2n,R) such that 

u(S) <Ip>* u(s)- 1 = s < ip >* . (i4) 


From linear algebra and group theory, we know that there is a unique polar de- 
composition S — RP for any element S in Sp(2n,R), where R is orthogonal, P is 
symmetric and positive definite, and both R and P are in Sp(2n,R) . Therefore we 
can always put S = exp(mp) exp(mp), where R = exp(mp), P = exp(mp), and 
both mjj and ms are elements in sp(2n,r) (mp is symmetric and unique, while 
mp is anti-symmetric and not unique) [11]. The element U(S ) E Mp(2n,R) 
which is unitary and satisfies ( 14 ) can be constructed as follows: 

U(S) = exp[4>(-mp)]exp[4(-mp)], ( 15 ) 
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where exp[<t>(— m P )] corresponds to a generalized squeezing and exp[<f>(— mp)] 
to a rotation in 2n-dimensional phase space. This decomposition is crucial in 
the construction of multimode SqCS since the Fock vacuum is an eigenstate of 
exp[ < l(— mp)], hence only the degrees of freedom in exp[$(— mp)] are effective in 
the SqCS constructed as C/ (-S') |0) [8]. 


Lemma [12] 

If M is a symmetric and positive definite 2n x 2n matrix, then there exists a 
matrix 5 € Sp(2n,R) (but not unique), such that 

"= s ‘(o !) 5 ’ (16) 


where uj = diag(u>i,cj 2 , . . . ,u;„), > 0 for all j. 

Remarks: 

(1) S G Sp(2n,R) if and only if S t JS = J by definition. 

(2) u>j is not an eigenvalue of M in general. 

(3) The eigenvalues of JM are itch’s, hence we can calculate w/s from JM 
as an ordinary eigenvalue problem. 

(4) If the matrix Cj corresponds to a 2— dimensional rotation on the 
< Xj y kj > plane, then 




(17) 


Therefore S in (16) can be replaced by CjS and hence is not unique. 


3 Unified Definition of Multimode Squeezed 
Coherent States 

The Hamiltonian we will use in this section is of n-mode, inhomogeneously 
quadratic, time-independent, and with its quadratic part positive definite: 

H = ^<lp>M<lp> t +V<H 5 >\ (18) 
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where M is a 2n x 2 n, symmetric and positive definite (hence invertible) real 
matrix and V is a 1 x 2 tx real vector. This kind of Hamiltonian can be transformed 
into the “standard form”: 


H = DoUqHoUq 1 D 0 1 + constant, (19) 

where D 0 is a displacement operator, U 0 is an operator in Mp(2n,R) and 

H 0 = J2ujiNi, wi> 0. (20) 


1=1 

Using the formulas discussed in last section, the derivation of (19) is straight- 
forward: 

H = \<U>M<U>' +V<U>* 

£ 

= ~[< V > + VM~ l ]M[< jp >‘ +M~ 1 V t ] + constant 
2 

= - < (|- x 0 ); (p - ko) > M < (q- x 0 )] (j> - ko) >* +constant 

2 

= -D 0 < p > Dq'MDo <Q',P> * A) -1 + constant 

2 

= A>[^ <q]p> M <q\p >‘]-Do 1 + constant 

= A,[J <q\P>S t \ n ” }S <q]P >*]i>o' 1 + constant 

2 V U U7 


, 1 „ . , 
= Do[~Uq < q\p > U 0 


= DoU 0 [^<q;p> 
A " 


u; 

. 0 tu 


:) 

° ^ U 0 < q\p >‘ U 0 ']D 0 1 + constant 
° ^ <q‘,P> t ]Ua l b o 1 + constant 


= A>#o[i £ ^( 9 ? + PiWo 'A, 1 + constant 

1 t=i 
n 

= + constant 


= D 0 UoH 0 U 0 1 D 0 1 + constant, 

where VM~ l = — < 3:0! >1 A) = D{< xq\ ko >) and Uo = U(S ). 


( 21 ) 
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Without loss of generality, we can always drop the constant term and consider 

H = DoUoHoUo'Do 1 - (22) 

It is easy to see that the normalized ground state of this Hamiltonian is: 

-CWoiO) oc Do exp[$(— mp)]|0), (23) 

which is a SqCS in general, it contains the coherent state ( Uq — 1, -Do 7^ 1) and 
the squeezed state (Uo ^ 1, Do = 1) as two special cases. 

Therefore we can take (23) as a unified definition of the multimode SqCS. 
However, since those u>;’s in Hq do not appear in (23), the correspondence between 
(22) and (23) is many-to-one. The non-uniqueness of S , hence Uo, will not cause 
any trouble, because we have shown that S is unique up to some 2-dimensional 
rotations in phase space, and rotations correspond to exp[$(— m#)] in U ( S ) which 
will not appear in (23). 

4 Multimode Thermal Squeezed Coherent States 

Consider immersing a physical system described by the Hamiltonian (22) in 
a heat bath of temperature T. This constitutes a canonical ensemble and the 
density operator of this system is: 

p = Z~ l exp(— PH) 

= Z-'expl-piDoUoHoUo 1 ^ 1 )] 

= Z~ l D 0 Uo ex-p(— (3 Ho)Uq 1 Dq 1 , (24) 

where 

/? = Z = Tr[exp(-/3H)] = Tr[exp(-/3H 0 )]. (25) 

This density operator p describes a mixed state unless T = 0. In the limit as 


T — » 0, since 

lim exp(-/3H 0 ) = |0)(0|, 

(26) 

we have 

p— >00 

^ = B o e'„|0){0|£) o - 1 C o -', 

(27) 
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which corresponds to the pure SqCS (23), hence (23) is a special case of (24) and 
(24) is a “thermalized state” of (23). Therefore we can take (24) as the definition 
of multimode thermal SqCS. 

5 Representations of Multimode Thermal Squeezed 
Coherent States 

There are many equivalent representations of the density operator p, e.g., 
the coordinate representation, P-representation, Q-representation, Fock space 
representation, Wigner function and characteristic function, etc. In this paper 
we will discuss the last two representations. 

5.1 Wigner Function 

The Wigner function of a density operator p is defined as [13, 14]: 

W(x\k) — 7 T~ n f dyexp(2ik • y)p(x — y,x + y), (28) 

J —oo 

where p(x, x') is the coordinate representation of the density operator p. 

The Wigner function can also be put into the following form [15]: 

W(x;£) = Tr[p£\ w {< x;k >)], (29) 

where the “Wigner operator” Aw(< x; k >) = 7r -n J9(2 < x; k >)exp(t7r^"_j TV,) 
is a well-defined hermitian operator with < x; k > as its parameters. 

The Wigner function is normalized by definition: 

/ oo _ 

dxdk\V(x;k) = 1, (30) 

-OO 

and it is real because the Wigner operator is hermitian. However, the Wigner 
function is not always positive-definite and it is thus called the quasi-probability 
distribution function over the “phase space” (x;£). 

In the following, we will calculate the Wigner functions for some thermal 
SqCS’s. First let us consider the simplest one-mode case, i.e., 

1 



the density operator is: 


(32) 


p — Z 1 exp(— /3u >N), 

and the Wigner function takes the form [14]: 

W(x, k) = Tr[pA w (x, *)] = ^ tanh(^) exp[- tanh(^y )(x 2 + k 2 )]. (33) 

In the limit as T —> 0, (33) becomes 

W(x, k) = — exp[— (x 2 + fc 2 )], (34) 

7 r 

which is exactly the Wigner function of the vacuum state [14]. 

Noticing that the Wigner function (33) is a Gaussian distribution function in 
(x, k ), we can use the exponent of (33) to define the “Wigner ellipse” in the phase 
space (x, k) as: 

tanh(^)(x 2 + k 2 ) = 1. (35) 

The area of the Wigner ellipse represents the range of uncertainty of the corre- 

sponding state. In this simplest case, the Wigner ellipse is a circle with radius 
^/coth(^/3o;) > 1 and with its center at the origin. 

Next we consider the general one-mode Hamiltonian: 

H = ujDoUqNUq 1 Dq 1 , (36) 

the density operator is: 

p = Z-'DoUoexpi-PuNWo'D^ 1 . (37) 


The Wigner function takes the form [16]: 


W(x,k) 


where 


Z l Tr[DoU 0 exp(—fiujN)U 0 1 D 0 1 Aiv(x, k)] 
Z~ l Tr[exp(— (3uN)Uq 1 Dq 1 Aw(x, k)D 0 Uo ] 


Z 1 Tr[exp(— /3uN)Aw(x' , fc')] 

-tanh(^)exp[-tanh(^)(x' 2 + k' 2 )], 
7T Z 4 



= S 


X - x 0 \ 

k - h ) ’ 


(38) 

(39) 
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and the Wigner ellipse is: 

tanh(^)(i — x 0 , k — koj&S ^ ^ = 1. (40) 

For the n-mode cases we first consider the uncoupled Hamiltonian, i.e., H = 
H 0 . The Wigner function in this case is a product of each individual one-mode 
Wigner function: 

W(x;k) = 7r" n [f[tanh(^ i )]exp[-X3 tanh (^ 1 )(^ + *?)]• ( 41 ) 

i=i z «=i z 

In this multimode case, we can define the “Wigner ellipsoid” in 2n-dimensional 
phase space as: 

f;tanh(^)(x? + *i) = l ) (42) 

1=1 ^ 

or equivalently, 

(^ ; ^) (q = ( 43 ) 

where T = diag(tanh( 4 /?wi),tanh( 4 /?a> 2 ), • • - ,tanh(|/?u> n )). Analogously, the 2n- 
dimensional volume of the Wigner ellipsoid represents the range of uncertainty. 

A A A A A _ A _ 

For the most general Hamiltonian H = DqUqHqUq Dq , analogue to (38), 
the Wigner function is: 

W(x]k) = 7r- n [JJtanh(^ 1 )] 

t=l Z 

exp[-(x-f 0 rf-k)? (J Q S(x-x 0 -,k-ko)% (44) 

and the Wigner ellipsoid becomes: 

(x — x 0 ]k — k)S* ^ ^ S(x — x 0 ]k — koY = 1. (45) 

5.2 Characteristic Function 

The characteristic function of a density operator p is defined as: 

X(x-k)=Tr[pD(-x--k)], (46) 
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From the symplectic Fourier transformation of the Wigner operator Au'(x; k) : 
F[Aw(x; k)] = [ dx'dk’Kw{x'\k')eyi^[-i{x' • k - k' • x)] 

= D{-x-,-k ), (47) 

we can see that the characteristic function is the symplectic Fourier transforma- 
tion of the Wigner function: 

F [W{x-k)} = X{x-k). (48) 

The normalization condition of the Wigner function corresponds to X(0; 0) = 
1 in the characteristic function. Since the operator D(—x] —k) is unitary instead 
of hermitian, X(x\ k ) is complex in general. 

The characteristic function of the general n-mode thermal SqCS, which cor- 
responds to the Wigner function (44), is: 



S{x-,kY + i{x-ko-k-x 0 )]. 


(49) 


5.3 Covariance Matrix 


For an n-mode (mixed) state w T ith density operator p, the covariance matrix 


is a 2n x 2n matrix of the form: 

V Q\ 

Q l v)' 

Uii = ((« - («»(& - («,-»> = (®9i) - <&>(&>» 

Vij s ((Pt - (Pi))(Pi - <Pi») = <PiPi> - 

Qn s ^ ((<?.- - <<?i>)(Pi “ (PiX + (Pj - <Pj ))(<?. - <&))) 

= ( \ (®Pi + Pj<?<)> - (<?*) (p? ) » 


(50) 

(51) 

(52) 


(53) 


where (<J,) = Tr(pqi), etc. 

For the thermal SqCS which corresponds to the Wigner function (44) or the 
characteristic function (49), it can be proved that the covariance matrix is 


1 

2 


S -1 


T -i 


0 


0 

T -1 


(S- 1 )*. 


(54) 
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6 Decoherence Problems 


6.1 General Theory of Decoherence 

Consider a quantum system which contains two subsystems (A) and (B) with 
the density operator Pab ■ Any (monomial) operator 0 which corresponds to a 
measurement on the system can be decomposed into Oa ® 0 B , where Oa corre- 
sponds to a measurement on and only on (A) and Ob correspondingly on (B). If 
we decohere this system by ignoring (B), i.e., not making any measurement on 
(B), then the operator 0 will be reduced to Oa <S> i and the expectation value of 
Oa will become: 


{Oa) = Tr[p A B(0 A ® i)] 

= T r {A )T r {B ) [ pab (O a <8> 1)] 

= Tr {A) [(Tr ( B)(pAB))d A ] 

= Tr[p A 0 A ], (55) 


where T r^ A )/ T' r (B) represents the “partial trace” which only takes trace with 
respect to the degrees of freedom of (A)/(B), and p A = Tr^ B ){pAB) is a well- 
defined reduced density operator. 

If the Wigner function W{x A , x B \ k A , k B ) corresponds to the original density 
operator Pab, then the reduced Wigner function corresponding to p A is [14]: 


_ r°° - _ 

W A {xA-,k A ) = / dx B dk B W{xA,x B -,k A ,k B ). 

J — OO 


(56) 


As for the characteristic function, if X{xa,x b \ k A ,k B ) corresponds to Pab, 
the reduced characteristic function corresponding to pa will take the form: 


X A {xA-,k A ) = X{xA,fr,k A ,Q), (57) 

which is a restriction of the original X{xA',x B ,kA‘,k B ) to a subspace in the 2 n- 
dimensional phase space. From the mathematical point of view, it is easier to use 
the characteristic function to study decoherence problems. 
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6.2 Decoherence of a Thermal Squeezed Coherent State 
From n — Mode to m— Mode 

For a given characteristic function of an n-mode thermal SqCS: 

X (x^ , X 2 , • • • ) X m , • • • , X n , k\ , k 2 , « • * } k m j ■ ' ' J ^n) 

= exp{— -(xi, x 2 , . . . ,x m , . . . ,x n ; k 1} k 2 , . . . , ■ ■ • ,^n) 


->-i 


0 T 


2=1 

The reduced characteristic function is: 


p_j ^ 5"(Xx j X 2 j • • • j X m , . . . , X n , k \ , fc 2 , • • • j • i ^-n) 


(58) 


^ (xi , x 2 > - ■ • ) x m , 0, /tj ? k 2 , • • • , Yn , 0) 

= exp{-^(x!,x 2 , . . . , x m , 0; fci, fc 2 , . . . , fc m ,0) 


5* 


---i 


0 


0 T 


_ j J S (x x , x 2 , . . . , x m , 0 , k \ , fc 2 , • • • •> k m , 0) 


^ (xjkpj kjXoj)} 


j=i 


— exp{ ^(Xj,X 2 , . . . , X m 5 fcj, fc 2 , • • • 5 &rn i 1 ) > ■ ■ ■ j ) ^-1 j ^-2 » • • • j ) 


+i'52(xihn - kiXoi)}, 

i=i 


(59) 


where the matrix I\ is 2m x 2m and still symmetric and positive-definite, its 

/T" 1 0 \ 

elements are a subset of the the elements of S l I _ x J S: 


( T~ l 0 \ 

•h i+m,j ~ [•S" ^ q j,_ x l 

/v, J+m = [5' (Y Y, 5 k , + „, 


(60) 

(61) 

(62) 
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(63) 


K, 





i+nj+n? 


where 1 < i, j, < in. 

From the Lemma in Sec. 
that: 


1 , we can find a 2m x 2m symplectic matrix <7 such 


I< = a 1 



(64) 


where r = diag(ri, 12 , . . . , t to ), t; > 0, for all i — 1,2, ... ,m. 

We can make a further restriction on t , from the following physical considera- 
tion: Since the reduced density operator p A = TV( B )(/>ab) is well-defined, it will 
never correspond to any non— physical state. Noticing that (59) is of the same 
form as (49), comparison with (54) shows the covariance matrix of this decohered 
state to be: 

Since c ~ * corresponds to a symplectic (hence canonical) transformation on the 
canonical coordinates, 



is also a covariance matrix for the same state in another canonical coordinates. 
This guarantees that r, > 1 for all i = 1 , 2 ,... ,m, otherwise (59) will give a state 
that violates the uncertainty principle. Therefore we conclude that the reduced 
characteristic function (59) corresponds to an m— mode thermal SqCS. 


0 




(65) 


7 Conclusion 

The results of this paper are threefold (1) A unified construction of multimode 
(thermal) SqCS’s. (2) Proof of the statement: The decohered multimode thermal 
SqCS is still a (multimode) thermal SqCS. (3) Introduction of the decohering 
technique via characteristic function, which is very efficient and can be applied 
to many related problems. 
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